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ABSTRACT 

We report on a measurement of the angular spectrum of the anisotropy of the microwave sky at 30 
and 40 GHz between I = 50 and I = 200. The data, covering roughly 600 deg 2 , support a rise in the 
angular spectrum to a maximum with STi w 85 /^K at I = 200. We also give a 2-sigma upper limit 
of STi < 122 /iK at I — 432 at 144 GHz. These results come from the first campaign of the Mobile 
Anisotropy Telescope (MAT) on Cerro Toco, Chile. To assist in assessing the site, we present plots of 
the fluctuations in atmospheric emission at 30 and 144 GHz. 

Subject headings: cosmic microwave background — cosmology: observations — atmospheric effects 



1. INTRODUCTION 

The characterization of the CMB anisotropy is essential 
for understanding the process of cosmic structure forma- 

, tion (e.g. Hu et al. 1997). If some of the currently pop- 
ular models prove correct, the anisotropy may be used to 

' strongly constrain cosmological parameters (Jungman et 
al. 1995, Bond et al. 1998). Here we report the results 

■ from the TOC097 campaign of the Mobile Anisotropy 
Telescope (MAT) experiment. 

2. INSTRUMENT 

The MAT telescope is comprised of the QMAP balloon 
gondola and instrument (Devlin et al. 1998), mounted 
on the azimuthal bearing of a surplus Nike Ajax military 
radar trailer 3 . The receiver has five cooled corrugated feed 
' horns, one at K a band (31 GHz), two at Q band (42 GHz), 
and two at D band (144 GHz). Each of the K a and Q band 

■ horns feed two HEMT-based (high electron mobility tran- 
sistor) amplifiers (Pospieszalski 1992) with one in each po- 
larization. The two D band horns each feed a single SIS 
detector (Kerr et al. 1993) with one horn in each polar- 
ization. This gives a total of eight radiometry channels 
in the experiment 4 . A Sumitomo mechanical refrigerator 
cools the HEMT amplifiers to 35 K and the SIS receivers 
to 4 K. 

The telescope optics are similar to those used for 
three ground-based observing campaigns in Saskatoon, SK 
(Wollack et al. 1997, SK). The feeds underilluminate 
an ambient temperature 0.85 m off-axis parabolic reflec- 
tor which in turn undcrilluminates a computer controlled 
1.8 mx 1.2m resonant chopping flat mirror. The beams 
are scanned horizontally across the sky in a ~s 4.6 Hz sinu- 
soidal pattern. The outputs of the detectors are AC cou- 
pled at 0.15 Hz and sampled N c times during each chopper 
cycle (N c — 80 for K a and Q bands, and N c = 320 for 
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2 University of Pennsylvania, Department of Physics and Astronomy, 



D band). The telescope is inside an aluminum ground 
screen which is fixed with respect to the receiver and 
parabola. 

The telescope pointing (Table 1) is established through 
observations of Jupiter and is monitored with two redun- 
dant encoders on both the azimuth bearing and on the 
chopper. The absolute errors in azimuth and elevation 
are 0?04, and the relative errors are <0?01. The chopper 
position is sampled 80 times per chop. When its rms posi- 
tion over one cycle deviates by more than 0?015 from the 
average position (due to wind, etc.), we reject the data. 



Table 1. TOC097 beam characteristics 



Feed 


Az 
deg 


El 

deg 


Qmeas 

10~ 4 sr 


nFWHM 
D avg 

deg 


K a l/2 


203.13 


41.75 


2.75 


0.90 


Ql/2 


206.75 


41.85 


1.69 


0.70 


Q3/4 


206.70 


39.25 


1.77 


0.72 


Dl 


205.00 


40.44 


0.183 


0.23 



3. OBSERVATIONS AND CALIBRATION 

Data were taken at a 5200 m site 5 on the side of Cerro 
Toco (lat. = -22?95 long. = 67?775 ), near San Pedro de 
Atacama, Chile, from Oct. 20, 1997 to Dec 15, 1997. The 
receiver was operational 90% of the available time. For 
the anisotropy data, the primary optical axis is fixed at 
az = 204?9, el = 40?5, 5 = -62?6 and the chopper scans 
with an azimuthal amplitude of 2? 96 (8? 93 on the sky) as 
the sky rotates through the beam. The telescope position 
was not wobbled to the other side of the South Celestial 
Pole as for the SK measurements in the North. The rms 
outputs of the K a 2 and Dl channels are shown in Fig. 1. 

NJ 08544 

David Rittenhouse Laboratory, Philadelphia, PA 19104 

3 Details of the experiment, synthesis vectors, data, and analysis code may be found at http://www.hep.upenn.edu/CBR/ and 
http:/ /pupgg. princeton.edu/~cmb 

4 HEMT amplifiers have improved considerably since this time (Pospieszalski et al. 1997) and SIS receivers are generally more sensitive than 
what we achieved. In 1997, one of the D channels and one of the Q channels did not have sufficient sensitivity to warrant a full analysis. 

5 The Cerro Toco site of the Universidad Catolica de Chile was made available through the generosity of Prof. Hernan Quintana, Dept. of 
Astronomy and Astrophysics. It is near the proposed MMA site. 

1 



2 



TOC097 



80 
70 
£ 60 

50 
40 

80 



50 
40 

















. ; 








~ w v : w' V" ^ \r' v >„' ' 












. *'- : i • ■ . ; c < * '• ■ : ■ : fc ; ' , ' , 
i'.yy y ^ Utf^y ^ ^ ^ H' ' V ^ i \t >' 


V* * V \i v 





290 300 310 320 330 340 350 



Fig. 1. — The rms detector output in antenna temperature of the K a 2 and DI channels averaged over 0.68 ms (with the chopper running 
at the nominal amplitude) vs. day of year in 1997. The sky is most stable between 10 PM and 10 AM local time. Similar results from 1998 
are consistent with the NRAO opacity measurements (http://www.tuc.nrao.edu/mma/sites/sites.html). 



Jupiter is used to calibrate all channels and map the 
beams. Its brightness temperature is 152, 160, 170 K for 
K a through D bands respectively (Griffin et al. 1986, 
Ulich et al. 1981), with an intrinsic calibration error of 
« 5%. We account for the variation in angular diameter. 
We also observe Jupiter with multiple relative azimuthal 
offsets to verify the chopper calibration. 

The uncertainty in the beam solid angle for the K a and 
Q bands is w 5% as determined from the standard devia- 
tion of beam measurements for the MAT and QMAP ex- 
periments. From a global fit of the clear-weather Jupiter 
calibrations, the standard deviation in the fitted ampli- 
tudes is 6%. These sources of calibration error dominate 
the error from the uncertainty in the passband. The total 
lcr calibration error is obtained by combining the intrinsic, 
beam, and measurement errors in quadrature resulting in 
10%, 10%, and 11% in K a through D respectively. 

A thermally-stabilized noise source at T c g « 1 K is 
switched on twice for 40 msec every 100 seconds as a rel- 
ative calibration. The pulse height is correlated to the 
Jupiter calibrations in the K a and Q channels. The vari- 
ation in detector gain corrected for with these calibration 
pulses is roughly 5%. No such correction was made for 
D band. 

4. DATA REDUCTION 

The reduction is similar to that of the SK experiment 
(Netterfield et al. 1997). The raw data, di, are multiplied 
by "n-pt" synthesis vectors, SV n j (where i ranges from 
1 to N c ) to yield the effective temperature correspond- 
ing to a multilobed beam on the sky, H(il). For exam- 
ple, we refer to the classic three-lobed beam produced by 
a "double difference" as the "3-pt harmonic" and write 
*3 = Efci SVz.idi- We also generate the quadrature sig- 
nal q n (data with chopper sweeping in one direction minus 
that with the sweeping in the other direction) and fast- 
dither signal /„ (one value of t n minus the subsequent 
one). For both K a 1/2 and Q3/4 we analyze the unpolar- 
ized weighted mean of the combined detector outputs. 

The phase of the data relative to the beam position is de- 
termined with both Jupiter and observations of the galaxy. 
We know we are properly phased when the quadrature sig- 
nal from the galaxy is zero for all harmonics. 

The harmonics are binned according to the right ascen- 
sion at the center of the chopper sweep. The number of 
bins depends on the band and harmonic (Table 2). For 



each night, we compute the mean and variance of all the 
t n , q n , and corresponding to a bin. These numbers are 
appropriately averaged over the campaign and used in the 
likelihood analysis. 

From the raw dataset of 814250 5s averages, we filter 
out time spent on instrument calibration (6%), celestial 
calibrations (11%), observations of the galaxy & daytime 
(53%), and bad pointing (4%). Accounting for overlap, 
these cut a total of 57%. The data span RA = 0° to 140° 
(b= -55° to -10°). 

The data are selected according to the weather by 
examining each harmonic independently. We first flag 
5s averages with a large rms. The unflagged data are di- 
vided up into 15 minute sections and the rms of the t n 
found. For 15 m sections with rms > 2a, the constituent 
5s averages are not used, as well as those of the preceding 
and succeeding 15 m sections. We ensure that the cut does 
not bias the statistical weight. As a final cut, nights with 
less than 4.7 hours of data are excluded. Repeating the 
analysis with increased cut values produces statistically 
similar (within lcr) results. The atmosphere cut selects 
roughly the same sections for K a and Q. In the analyses, 
we discard the 2 and 3-pt data as it is corrupted by atmo- 
spheric fluctuations and variable instrumental offsets. If 
the 4-pt is corrupted, it is at the la level and not readily 
detectable. 

The stability of the instrument is assessed through in- 
ternal consistency checks and with the distribution of the 
offset of each harmonic. The offset is the average of a night 
of data after the cuts have been applied (ranges from 5- 
10 hours) and is of magnitude w 200 /iK with error 20 /iK. 
In general, the offset remains constant for a few nights and 
then jumps 3-5 sigma. The resulting x 2 /v is typically be- 
tween 4 and 20 for the data over the full campaign and 
is w 1 for the quadrature signal. In general, a change in 
offset can have any time scale. The q n and are sensitive 
to r = 0.25 s. We also monitor a slow dither (difference of 
the subsequent 5 sec averages) with r = 5 s and a night- 
to-night dither with r = 24 h. For the final analysis, we 
delete one seven day section that has a large jump in off- 
set. To eliminate the potential effect of slow variations in 
offset, we remove the slope and mean for each night. This 
is accounted for in the quoted result (both the constraint 
matrix method, Bond et al. 1998b, and marginalization, 
Bond et al. 1991 give similar corrections) and does not sig- 
nificantly alter the results over the subtraction of a simple 
mean. As a test, we have also tried removing quadratic 



Torbet et al. 



3 



and cubic terms from the offset, with no significant changes 
in the answer. In summary, there is no evidence that the 
small instability in the offset affects our results. 

We examined the variations in the power spectrum of 
the synchronously co-added raw HEMT data, and found 
no evidence for microphonics. However, a microphonic 
coupling to the SIS detector was exacerbated after situ- 
ating the telescope at the site. After filtering, residual 
signals persisted in the quadrature channels (though not 
in the fast and slow dithers) and so we report only 95% 
upper limits for the D channel, specifically 6Ti < 180 fiK 
at£ = 325 and 5Ti < 122 at i = 432. 

The primary effect of data editing is to increase the er- 
ror bar per point and decrease the upper limits of the null 
tests. Of the 169 null tests (Table 2 plus fast, slow, and 
night dithers), there are only three failures. The distribu- 
tion of the reduced y 2 of the null tests is consistent with 
noise and inconsistent with any signal. When the data 
are combined into groups of harmonics and bands, all null 
tests are consistent with noise. 

5. ANALYSIS AND DISCUSSION 

The analysis of the individual harmonics, because the 
windows are so narrow, essentially corresponds to finding 
6T{ = v/(A| ot - A? nat )/J(W0, where A t 2 ot is the variance 



of the covariancc matrix. 



of the data for each harmonic, A 2 nst is the variance due to 
atmospheric and instrumental noise, and I(W) = Wi/l. 
Wi is the window function, as defined in Bond 1996. The 
full likelihood analysis provides a formal way of determin- 
ing STi that includes correlations and gives the correct 
error bar in the low signal-to-noise limit. 

The error in I(W) is determined from the scatter in the 
beam values. We find 6I(W)/I(W) £ 0.01 for all bands 
and harmonics. The mean variance, A 2 nst , is determined 
directly from the uncertainties in each bin. If these un- 
certainties are somehow biased, the results of the simple 
test and full likelihood will be biased. We examine the 
distribution of all the data for each harmonic from all the 
nights after removing the mean value of each sky bin. The 
width of this distribution agrees with the mean error bar 
indicating that the error per point is not biased. Also, the 
ratio of the error bars between harmonics agrees with the 
analytic calculation. 

In the full analysis (Fig. 2), we include all known corre- 
lations inherent in the observing strategy. From the data, 
we determine the correlations between harmonics due to 
the atmosphere, detector noise, and non-orthogonality of 
the synthesis vectors. The correlation coefficients between 
bands due to the atmosphere are of order 0.05. We also 
examine the autocorrelation function of the data for a sin- 
gle harmonic to ensure that atmospheric fluctuations do 
not correlate one bin to the next. The quoted results are 
insensitive to the precise values of the off-diagonal terms 



ll 



\\j 



Fig. 2.- 

values are 
,1+20 V n+6 



Combined analysis of data in Table 2 
(l,6Ti [/iK]) = (63+ 18 " n+10 ^ ^«+ 16 



-18i w -9 



The 

(86±||,45±S), 



(114+^,70^), (158±^,89±£), (199±i*,85tf). Error bars 
are "lcr statistical"; calibration error is not included. The 
COBE/DMR points are from Tegmark 1997. The solid curve 
is standard CDM (Q h = 0.05, h = 0.5). 

These results are similar to previous results obtained 
with this technique (SK) though the experiment was done 
with different optics, a different receiver, a different pri- 
mary calibrator, largely different analysis code, and ob- 
served a different part of the sky. Though we have not 
correlated our data with templates of foreground emission, 
the foreground contribution is known to be small at these 
frequencies and galactic latitudes (Coble et al. 1999, de 
Oliveira-Costa et al. 1997). In addition we have exam- 
ined the frequency spectrum of the fluctuations in K a and 
Q bands, and find it to be consistent with a thermal CMB 
spectrum, and inconsistent with various foregrounds. Fi- 
nally, the full analysis has been repeated after deleting 
each 15° section of data in RA, indicating that the signal 
does not arise from one region. (Our scan passes near, but 
misses, the LMC.) Future work will address the precise 
level of contamination. 
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Table 2 
TOC097 Angular Spectrum 



Band/SV 


^cff a 


ST e b 


-^V bins 


ST' e A to t 






(A-B)/2 d > c 


Quad, q„ d 


(Hl-H2)/2 d - f 












fiK 


//K 







K a l/2 



4pt 


63±1I 


35+ 9 3 


48(16) 


32 


33 


20 


0.84 


< 27(0.94) 


< 30(1.05) 


< 29(0.95) 


5pt 


86i« 


52+j 1 


64(28) 


49 


40 


21 


0.71 


< 21(0.59) 


< 32(1.23) 


< 29(1.05) 


6pt 


107t« 


n +12 
"-lO 


96(42) 


69 


52 


27 


0.65 


< 32(0.90) 


< 31(0.96) 


< 35(1.04) 


7pt 




93+_\l 


96(41) 


90 


57 


27 


0.55 


< 35(0.80) 


< 37(0.90) 


< 30(0.68) 


8pt 


145+18 


103111 


128(55) 


102 


63 


34 


0.52 


< 51(0.97) 


< 46(1.00) 


< 43(0.91) 


9pt 


i65+i 8 




128(54) 


59 


47 


38 


0.46 


< 63(0.95) 


< 72(1.19) 


< 66(1.07) 


lOpt 


182+ji 




192(82) 


70 


60 


51 


0.44 


< 68(0.85) 


< 69(0.94) 


< 60(0.89) 


llpt 


192+f 


< 119 


192(82) 


67 


65 


58 


0.42 


< 91(0.94) 


< 83(1.00) 


< 86(0.96) 


12pt 


215+H 


128l|§ 


192(82) 


127 


83 


68 


0.37 


< 150(1.11) 


< 86(0.79) 


< 76(0.67) 



Qi 



4pt 


63±H 


57 +i8 

°'-13 


48(20) 


51 


53 


31 


0.83 


< 44(1.04) 


< 47(1.11) 


5pt 


87l 2 6 2 


40li 4 4 


64(28) 


34 


40 


33 


0.71 


< 36(0.75) 


< 47(1.09) 


6pt 


HOl 24 




96(42) 


52 


52 


40 


0.65 


< 45(0.87) 


52115(1.71) 


7pt 




01 +19 
81-16 


96(42) 


77 


59 


41 


0.55 


< 65(1.12) 


< 53(0.89) 


8pt 


15lt 2 5 5 


8611? 


128(55) 


79 


66 


50 


0.53 


< 45(0.67) 


< 64(0.89) 


9pt 


172+2* 


qo+23 


128(55) 


89 


68 


54 


0.47 


< 82(0.94) 


< 72(0.83) 


lOpt 


1Q1+16 
lyl - 2 4 


< 115 


192(84) 


31 


75 


74 


0.47 


< 105(1.13) 


< 92(0.94) 


llpt 


203l 2 7 


< 117 


192(84) 


44 


80 


78 


0.44 


< 103(0.94) 


< 83(0.81) 


12pt 


221 +2 j 


< 169 


192(84) 


91 


91 


84 


0.40 


< 138(1.06) 


< 119(0.92) 


13pt 


245l 2 2 J 


< 130 


192(84) 




84 


92 


0.37 


< 163(1.02) 


< 164(1.06) 


14pt 


267+_f 3 


< 202 


256(112) 


123 


123 


114 


0.37 


< 183(1.06) 


< 188(1.00) 



Q3/4 



5pt 


go+15 
83_ 20 


47+17 


64(17) 


43 


39 


24 


0.73 


3llg°(1.94) 


< 50(1.39) 


< 58(1.73) 


6pt 


106l 2 4 


61 +18 


96(22) 


55 


46 


28 


0.66 


< 27(0.68) 


< 34(0.60) 


< 61(1.46) 


7pt 


125l 2 4 


79 + 16 
' Z -12 


96(35) 


71 


50 


30 


0.56 


< 29(0.65) 


< 57(1.30) 


< 31(0.57) 


8pt 


145l 2 4 


11511? 


128(45) 


109 


68 


35 


0.54 


< 26(0.50) 


< 38(0.73) 


< 60(1.18) 


9pt 


165l 2 4 4 


79+24 
,z -21 


128(29) 


65 


48 


37 


0.48 


< 43(0.82) 


< 97(1.38) 


< 84(1.16) 


lOpt 


184l 2 5 3 


8711? 


192(70) 


86 


63 


48 


0.47 


< 68(1.03) 


< 45(0.78) 


< 52(0.75) 


llpt 


1961 22 . 




192(54) 


84 


65 


53 


0.44 


< 64(0.90) 


< 78(0.91) 


< 75(0.78) 


12pt 


212+H 


tool 2 ? 


192(65) 


97 


69 


57 


0.40 


< 70(0.95) 


< 129(1.36) 


< 84(0.89) 


13pt 


2361 fa 


< 157 


192(56) 


80 


71 


65 


0.37 


< 107(1.07) 


< 121(0.93) 


< 153(1.17) 


14pt 


258l 2 * 


1 1 Q+ 36 


256(103) 


119 


91 


80 


0.36 


< 116(0.95) 


< 137(1.11) 


< 103(0.78) 



Note. — A "<" indicates a 95% confidence limit. Calibration errors are not included, (a) The range for <? e g denotes the range for 
which the window function exceeds e~l/ 2 times the peak value, (b) The error on STp = [1(1 + l)C^/27r]l/ 2 is comprised of experimental 
uncertainty and sample variance. These values are not statistically independent: harmonic numbers differing by 2 are correlated at the 
0.35 level. For all harmonics, the sample variance (oc l/\/2Ni >i „ a ) is ss 7 fiK. (c) The number of bins on the sky followed by, in parentheses, 
the number used in the analysis due to the galactic/atmosphere cut. (d) The reduced x 2 are given in parentheses, (c) (A — B)/2 is the 
difference in polarizations. We have combined bands and harmonics to generate 95% upper limits on polarization, A — B, and obtain: 
(/, STi luK]) = (6311*, < 3 7): (86l 2 l, < 54), (1151^, < 28), (1481^, < 41), (195l||, < 79). (f) (HI - H2)/2 is the first half minus the 
second half. 
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Fig. 3.— Angular spectrum from the SK, QMAP, and TOC097 experiments. The SK data from Netterfield et al. (1997, ApJ, 474:47) and 
Wollack et al. (1997, ApJ 476:440) have been recalibrated according to the Mason et al. (1997, astro-ph/9903383), leading to an increase of 
5%, and reduced according to the forground contribution in de Oliveira-Costa et al. (1997, ApJL, 482:L17), leading to a reduction of 2%. The 
foreground reduction was applied uniformly. Because the SK data are primarily at 40 GHz, the Mason et al. results must be extrapolated. In 
addition, uncertainty in the measured passband, measurement uncertainty, and beam uncertainty, must be taken into accounted. This leads 
to a lcr calibration error of 11%. The QMAP data are the same as those reported in de Oliveira-Costa et al. (1998, ApJL, 509:L78) and have 
an average calibration error of 12%, the correction for foreground emission is PS 2% (de Oliveira-Costa et al. 1999, in prep), though has not 
yet been precisely determined and so is not included. Both SK and QMAP are calibrated with respect to Cas-A. The TOC097 data, which 
have a calibration error of 10%, are calibrated with respect to Jupiter. A foreground contribution, which is expected to be small, has not been 
subtracted. When calibration errors are included, all three independent experiments agree. Two cosmological models are shown for reference. 
The lower is "standard CDM" (f2o = 1, !2b = 0.05, h = 0.5); the higher one is the current "concordance model" (Bachall et al. 1999 in prep, 
Turner 1999 astro-ph/9904051) with Q = 0.4, fl B = 0.045, O a = 0.6, and h = 0.65. For COBE/DMR we use the results from Tegmark. 



